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ABSTRACT 

We report the first direct observations of neutral, molecular gas streaming 
in the nucleus of NGC1068 on scales of < 30 pc using SINFONI near-infrared 
integral field spectroscopy. At a resolution of 0.075", the flux map of 2.12 /iin 
1-0 S(l) molecular hydrogen emission around the nucleus in the central arcsec 
reveals two prominent linear structures leading to the AGN from the north and 
south. The kinematics of the gas in these features are dominated by non-circular 
motions and indicate that material is streaming towards the nucleus on highly 
elliptical or parabolic trajectories whose orientations are compatible with that of 
the disk plane of the galaxy. We interpret the data as evidence for fueling of gas 
to the central region. The radial transport rate from ~ 30 pc to a few parsec 
from the nucleus is ~15 M yr" 1 . One of the infalling clouds lies directly in front 
of the central engine. We interpret it as a tidally disrupted streamer that forms 
the optically thick outerpart of an amorphous clumpy molecular /dusty structure 
which contributes to the nuclear obscuration. 

Subject headings: galaxies: active - galaxies: individual (NGC 1068) - galaxies: 
nuclei - galaxies: Seyfert - galaxies: molecular gas - infrared: galaxies 



1 Based on observations at the European Southern Observatory VLT (076.B-0098). 
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Introduction 



A detailed description of the distribution and kinematics of the molecular gas in the 
central region of Seyfert galaxies is crucial for understanding the fueling of the nucleus and 
the ro le of gas in obscuring the acti ve galactic nucleus (AGN). NGC1068, at a distance of 14.4 
Mpc (iBland-Hawthorn et al .11199 71 ) . is the brightest Seyfert 2 galaxy and therefore one of the 
best candidates for direct investigations of the morphology and dynamics of the molecular 
gas neighboring an active nucleus. The idea of a rotating molecular /dust y torus surrounding 



a sup ermassive black hole emerged from the interpretation given by iMiller fc Antonucci 



( 119831 ) of optical spectropolarimetry of precisely this galaxy, which revealed scattered type-I 
emission from the obscured broad-line region (BLR). Since then, several near- and mid- 
IR observations have in fact found and begun to elucidate the physical conditions of the 
concentration of molecular gas and dust in the nucleus of NGC1068 (near-IR: Young et al. 
1996; Marco et al. 1997; Rouan et al. 1998, 2004; Alloin et al. 2001; Galliano & Alloin 2002; 
Galliano et al. 2003; Weigelt et al. 2004; Wittkowski et al. 2004; Gratadour et al. 2005, 
2006; mid-IR: Bock et al. 2000, Tomono et al. 2001, 2006; Jaffe et al. 2004; Galliano et al. 
2005; Mason et al. 2006; Poncelet et al. 2006, 2007), all favouring indirectly the existence of 
a compact molecular /dusty torus but failing to obtain a cl ear image of it. Altern atively, by 
means of mid-IR observations over the central 140 pc (2"), ICameron et al.l (119931 ) proposed 
a model in which the bulk of the molecular gas and dust is located at large distances (several 
tens of parsecs) from the AGN. Line-of-sight attenuation of the BLR in this cas e would be 



a mer e consequence of one or more intervening molecular clouds. More recently, iJaffe et al. 



( 120071 ) analyzed new interferometric mid-IR observations of NGC1068 and found that the 
fitted Gaussian components to the (u, v) plane of the central 10 /im source resemble a disk 
similar to the H 2 masers (jGreenhill et al.lll996l ). Both, the dust and H 2 maser disks, 
appear to be oriented neither perpendicular to nor aligned with the radio jet. However, one 
needs to be cautious when interpreting all of these observations. As they are measurements 
of the continuum emission, the inferred gas/dust distributions are strongly dependent on 
temperature, tracing in fact radiation of matter at a given temperature rather than spatial 
gas/dust distributions. 

The 2.122 /mi FLj rovibrational v =1-0 S(l) emission line probes hot (> 10 3 K) and 
moderately dense (> 10 3 cm -3 ) molecular gas and as such, may be an excellent tracer of gas 
in the nuclear region. Its spatial distribution can expose the potential presence of a molec- 
ular/dusty torus, and also provide physical information on fueling or feed back mechanisms. 
Imaging spectroscopy of this line (IRotaciuc et al.lll99ll ; iBlietz et al.lll994j ) has indicated the 
presence of significant amounts of hot, dense, circumnuclear molecular gas extending over 
the central 4 arcseconds which is associated with the narrow-line clouds. The nuclear re- 
gion shows a strong peak almost 1" east of the nucleus and a weaker one ~ 1" to the west. 
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Furthermore, the H2 emiss ion at a few arcsec onds from the n ucleus is more exte n ded a long; 
the major axis of the bar (jDavies et al.lll998l ). More recently, iGalliano fc Alloinl (120021 ) ob- 
tained 2D spectroscopic observations of the warm molecular gas in the central 4" x 4" of 
the galaxy. They did not detect H 2 emission at the location of the K— band continuum core 
(which is known already to be coincident with the central engine at these scales). Instead 
they confirmed the two main regions of H2 emission at about 1" east and west of the nucleus 
along a PA=90°, and a region with complex line profiles at ~ 50 pc north of the AGN. They 
interpreted the observed H 2 em ission at these scales in terms of a warped disk consistent 
with the interpretation given by ISchinnerer et al.l (120001 ) to millimeter interferometric maps 
of 12 CO(2-l) emission. It is now apparent from recent SINF ONI H 2 1-0 S(T) data at these 
scale s that simple warped disk models of the molecular gas (ISchinnerer et al.l |2000| ; [Baker 
20001 ) cannot account for the fantastic variety and detail in the morphological and kinemat- 
ical structure (see Davies et al. 2006), in particular since there is a considerable amount of 
gas inside the inner edge of the 100 pc ring. 

In order to analyze in detail the physical conditions of the warm molecular gas in the 
central region of an active nucleus, we have obtained adaptive optics assisted SINFONI 
integral field spectroscopic data to peer deeply into the energetic core of NGC1068. In this 
paper we focus on the observations of the molecular gas emission in the central arcsecond of 
the galaxy done with the smallest spatial scale of the instrument, resulting in a resolution 
of 75 milli-arcseconds, approxim ately 7 times better than previous observations. We refer 
to iMiiller Sanchez et al.l (l2008bl ) for the analysis of the stars, molecular and ionized gas in 



the central 3" of the galaxy corresponding to the rest of the results from the SINFONI 
observations of NGC1068. 



2. Observations 

The data presented here were obtained during 2005 -2006 using the adaptive optics as- 



sisted near-infrared integral field spectrograph SINFONI (IBonnet et al.ll2004l ; lEisenhauer et al 



20031 ) on the VLT UT4. SINFONI delivers spectra simultaneously over a contiguous two- 
dimensional field of 64 x 32 pixels. The data were taken in the H + K— bands using two of 
the three possible spatial pixel scales of the instrument: 0.0125" x 0.025" and 0.05" x 0.1". 

The 0.0125" x 0.025" SINFONI data were taken in two sets of 2 and 1.5h integrations 
on the nights of 21 October 2005 and 27 November 2006, both using the galaxy nucleus as 
reference for the MACAO adaptive optics system. In both nights the atmospheric conditions 
were excellent (optical seeing of the two nights was ~ 0.7" and ~ 0.5" respectively), which 
allows us to obtain a resolution of 0.075" FWHM as measured from the spatially marginally 
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resolved non-stellar continuum in if— band. The 0.05" x 0.1" data were obtained on 6 
October 2005. Once more, the AO module was able to correct on the nucleus of NGC1068 
in seeing of 0.63", to reach ~ 0.1" spatial resolution. In this total of 4 sky and 8 

on-source exposures of 50 sec each, were combined to make the final data cube. For the 
0.0125" x 0.025" and 0.05" x 0.1" datasets, the spectrograph was able, in a single shot, to 
obtain spectra covering the whole H and if— bands (approximately 1.45 — 2.45 /im) at a 
spectral resolution FWHM of 125 km s" 1 (R~2400) for each pixel in the 0.8" x 0.8" and 
3.2" x 3.2" field of view. 



The data were reduced using the SINFONI custom reduction package SPRED flAbuter et al 



20051 ). This performs all the usual steps needed to reduce near infrared spectra, but with 
the additional routines for reconstructing the data cube. Flux calibration was performed 
simultaneously in H and if— bands using a G2V star HD20339 (H = 6.45, if =6. 41), which 
yielded in a 1" aperture a if —band magnitude for NGC1068 of 7.73 and H— band magnitude 
of 10.18. The values we dedu ced are consistent, at similar and smaller apertures, with va l- 
ues found in previous studies (IRouan et al.lll998l ; iGratadour et al.ll2006l ; iPrieto et al.l 120051 ). 
In addition, flux calibration was cross-checked with VLT NACO data in 2 — 4" apertures 
( iMiiller Sanchez et al.ll2008bl ). The NACO data was further cross-checked in larger 5 — 10" 
apertures using 2MASS data. Agreement between the different data sources was consistent 
to 15%. 



3. Distribution, kinematics and physical properties of the molecular gas 

3.1. Molecular gas morphology 



The molecular hydrogen emission at scales of a few arcseconds from the nucleus has been 
mapped previously reaching spatial resolutions down to ~ 0.5" (Galliano & Alloin 2002). 
Our SINFONI data at these scales reach ~ 0.1" resolution and reveal a complex distribution 
of the gas which has not been entirely observed before (Fig. [1]). The new data resolve the 
previously studied H2 knots and show the presence of an off-center (Ar = 0.6" SW from the 
AGN) ring of molecular gas (see Fig. [H]), as well as apparently linear gas streamers leading 
from the ring NNW and SSE to the center . In ad dition, the overlay of our H 2 flux map and 
the 12 CO(2-l) map from ISchinnerer et al.l (120001 ) at these scales (Fig. [1]) shows a very good 
correlation, confirming that gas distributions can be traced either by the 2.122 fim H2 1-0 
S(l) emission line (hot gas) or the 12 CO(2-l) emi ssion (cold gas). We discuss briefly the 
morphology of the ring in Section 13.61 but refer to IMiiller Sanchez et al.l ( )2008bl ) for more 
analysis and discussion of these results. In the following we will concentrate on the central 
arcsecond (white box in Fig. [1]) and provide a few remarks on the connection of the linear 
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streamers to the circumnuclear environment. 



The morphology of the H 2 1-0S(1) emission in the central 0.8" x 0.8" of NGC1068 with 
a resolution of 0.075" (~5 pc) is presented in Fig. [2J The peak of the 2.1 /mi non-stellar 
continuum is located at the origin of the image and is represented b y a crossed circle. A t 
these scales this can be identified as the position of the central engine ([Galliano et al.ll2003l ). 
At larger scales, a linear structure leading to the AGN from the NNW and SSE is the only 
noticeable feature in this whole region (see Fig. [I]). On sub-arcsecond scales near the AGN, 
this linear structure exhibits two prominent regions of H 2 emission, one overlapping the 
center of the non-stellar continuum, and another one located ~0.35" north of this point. 

Thanks to improved spatial resolution and sensitivity, molec ular gas is now detecte d 
closer and closer to the AGN, in contrast to previous observations (jGalliano fc Alloinll2002l ). 
In Fig. [31 we show a spectrum of NGC1068 at a spectral resolution of R~2400, integrated 
over a 0.3" x 0.1" rectangular aperture centered at the position of the nucleus; although the 
continuum is quite strong in this region, the H2 1-0 S(l) line is clearly seen. We will refer to 
this nuclear concentration of gas as the southern tongue. It has a major axis of ~17 pc and 
a minor axis of ~7 pc and is elongated to the south following the SSE extension of the linear 
struc ture. Its major axis ha s a position angle of about 120°, similar to that of the H2O maser 



disk iGreenhill et al.l [l996l l and the 10 /im dust emission (jJaffe et al.l 12004 ; iPoncelet et al.l 



20071 ) on a scale of 20 mas. These three features exhibit the same misa llignment with the 
major axis of the nuclear 5 GHz radio continuum (iGallimore et al.l 12004 ) . The H 2 1-0 S(l) 
luminosity of the southern tongue is 1.5 x 10~ 18 W m -2 (~lx 1O 4 L ). 

In addition to the nuclear component designated as the southern tongue, we found 
another area of high luminosity in the northern part of the H2 intensity map (see Fig. [2]). 
In fact, this is the brightest region of the linear structure (~ 1.6 x 1O 4 L ) which extends 
up to ~ 0.7" NNW of the central engine along a PA of —14° and apparently connects to a 
prominent knot in the circumnuclear ring (Fig. [1]). We will designate this extension of the 
linear structure as the northern tongue. The region north of the AGN where apparently 
the nort hern tongue and the circ um-nuclear ring merge is characterised by double-peaked line 
profiles JGalliano fe Allomll2002h . As this region is out of the central arcsecond, we po stpone 
the quantitative analysis of it to a subsequent publication (iMiiller Sanchez et al.ll2008bl ). The 
H2 line profiles of the small scale data studied here do not show complex morphologies and 
thus they were fitted by one Gaussian component as described in Sect. 13.21 



Dust emission at 12 



um fr o m the norther n tong ue is present in pre v ious m id-IR ob 



servations by iBock et al.l (J2000j) , iTomono et al.l (120061 ) and IPoncelet et al.l (120071 ) . Overlays 
of our two H 2 flu x maps at different spatial scales and the 12.5 /im deconvolved image of 
Bock et al.l (120001 ) are shown in Fig. HI The image taken at 12.5 jum has the highest angular 
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resolution and highest signal to noise ratio of the mid-IR images from iBock et al.l (120001 ) . 
For these reasons we have chosen to compare this image with our flux maps, in a similar way 
as these authors did with other datasets. We want to emphasize, however, that all mid-IR 
images, show qualitatively the same features. In both panels of Fig. @] we assumed that 
the mid-IR peak is located at the position of the nucleus as defined by the near-IR peak. 
As can be seen in the two overlays, along the northern tongue there exists good correlation 
between the molecular gas emission and the 12.5 fim continuum. In the right panel, a dashed 
line indicates the boundary of the northern tongue emphazising this correlation. Also, the 
alignment gives a fair correlation between the east-west unresolved mid-IR core and the H2 
southern tongue. The two overlays show a bend to the east, although the mid-IR bend is 
more pronounced. The qualitative agreement of the two images indicates that in the central 
arcsecond the molecular gas and the dust have a similar spatial distribution which is pre- 
dominantly a north-south linear structure about 1" (70 pc) long which contains two bright 
components: the southern and nor thern tongues. Our observations support the interpreta- 



tion given by iTomono et al.l (120061 ) suggesting that the northern tongue could be a passage 
that transports material to the central engine. We show for the first time, direct evidence 
from the gas kinematics that this is indeed what is occurring (see Sect. 



Finally, a comparison of the H 2 and 5 GHz flux maps is shown in Fig. [5j Radio continuum 
imaging at 5 GHz with a resolution of 0.065" revealed a number of structures along the inner 
part of the radio jet (Gallimore et al. 1996). Associating component SI with the inner edge 
of the torus, these authors developed a scenario in which component C arises from a shock 
interaction between the jet and a dense molecular cloud. Supporting this hypothesis was the 
bend in the radio jet, the slightly flatter spectral index, and the presence of maser emission. 
Fig. [5] shows that component C does in fact coincide spatially with the northern tongue of H 2 



1-0 S(l) emission. Moreover, the u = 1 levels are thermalised (jMuller Sanchez et al.ll2008bl ). 
indicating that the gas is likely to be rather dense (> 10 4 cm -3 ). Thus the SINFONI data 
provide direct evidence for the molecular cloud and hence strongly support the jet-cloud 
interaction hypothesis. 



3.2. Kinematic evidence for inflow 

3.2.1. Detection of non- circular motions around the nucleus 

We have extracted gas kinematics in the nuclear region of the galaxy from our integral 
field data, allowing for the first time a study of the gas motions on scales connecting the 
outer gas r ~ 1" (~70 pc) ring to the maser disk at r = 0.015" (~1 pc). Fig. E] shows the 
velocity and dispersion maps of the molecular gas in the central 0.8" x 0.8" of the galaxy. 
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We have extracted the 2D kinematics by fitting a Gaussian convolved with a spectrally 
unresolved template profile (an OH sky emission line) to the continuum-subtracted H 2 1-0 
S(l) spectral profile at each spatial pixel in the data cube. We performed a minimisation 
of the reduced \ 2 m which the parameters of the Gaussian (amplitude, center and width 
in velocity space) were adjusted until the convolved profile best matched the data. The 
uncertainties were boot-strapped using Monte Carlo techniques, assumi ng that the noise i s 



uncorrelated and the intrinsic profile is well represented by a Gaussian (jDavies et al.l 120071 ). 
This method allowed us to obtain uncertainties for the velocity and dispersion in the range 
of ±(5 — 15) km s _1 . The dispersion extracted by this fitting procedure is already corrected 
for instrumental broadening. 

The velocity field shown in Fig. [H] is quite complex. There is no evidence of a rotating 
disk around the nucleus. Instead, one can distinguish three main kinematical components: 
blue-shifted material with almost constant projected velocity v z of —25 km s -1 in the north, 
red-shifted projected velocities between 20 — 40 km s -1 in the south-east, and a nuclear red- 
shifted component (v z ~ 90 km s~ x ) associated with the southern tongue in Fig. [2j Outside 
these regions the kinematics are not reliable due to the low strength of the line (< 10% of 
the flux as can be seen in Fig. [2]) and therefore they were masked out in the velocity and 
dispersion maps. The nuclear red-shifted kinematical component appears to be connected 
to the south-eastern kinematical component by a ridge of emission that gradually changes 
its projected velocity from ~ 30 km s _1 at r = 0.4" to ~ 90 km s _1 at r = 0.04" along 
a PA of —20°. This and the lack of a signature of rotation indicate that the gas must be 
streaming almost directly towards or outwards from the nucleus rather than orbiting it on 
circular paths. Thus, any approach to reproduce the observed kinematics by a rotating or 
warped disk model can be excluded. 



3.2.2. Quantitative modelling 



To quantify the kinematics of the region, we have modeled the observed velocities as 
motions of test particles under a gravitational potential comprising a central mass and an 
extended stellar component. 

We define a three-dimensional cartesian coordinate system in the nuclear region of the 
galaxy with x—, y— and z— axes representing the Right Ascencion, Declination and Line-of- 
Sight directions respectively. The gravitational potential well of the system is determined 
by a circum-nuclear mass distribution for med by the sum of a supermassive black hole of 
M BH = 1 x 10 7 M flGreenhill et al.lll996f ) located at th e origin of the system, and a stella r 
mass density M*(r) = 1 x 10 6 r M pc _1 JPavies et aDl2007l : Muller Sanchez et al.ll2008bh . 
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The gas clouds were modelled as test particles - i.e. they do not have any impact on the 
potential - with some initial conditions x , y , z , v x0 , v y0 and v z0 . The initial x , y and 
v g Q components of a test particle are obtained directly from the SINFONI data. Thus, 
the kinematic model contains three degrees of freedom corresponding to the other three 
Cartesian components z Q , v x q and v y o- Notice that v z q can be considered also a variable due 
to the uncertainity in the observed projected velocities of ±10 km s -1 . Once the phase space 
conditions of a cloud are given, the position and velocity of the cloud moving according to 
the Newtonian laws of motion under the influence of the assumed potential were determined 
at every time interval At = 10 yr over a period of 5 Myr. Therefore, this method creates free 
unperturbed Keplerian orbits fully defined by the initial values of the Cartesian components 
of the particle's position and velocity. 

We followed a systematic approach for the determination of the Keplerian orbits of the 
gas particles. First, an initial position vector ro is located inside a volume defined by the 
field of view and several Zq components ranging from —60 pc to 60 pc. This interval was 
defined as two times the x— or y— range. It is important to point out that this volume 
delineates the boundaries of the initial position vectors. The resulting orbits are actually 
contained in a spherical volume of indefinite radius. For each r numerous sets of velocity 
vectors were investigated. The initial v z q component of a particle at any given position 
in space corresponds basically to the observed projected velocity at that particular spot. 
The tested values for v x o and v y o at each point in space were dependent on r and ranged 
from — a/2GMbh/^o to a/2GMbh/?"o- This interval was established based on the universal 
expression for the tangential velocity v tan of any point on the orbit. If the magnitude of the 
initial tangential velocity vector is smaller than this factor, the motion will be elliptic, if it 
is larger, the motion will be hyperbolic, and if it is precisely this value, a parabolic orbit 
will be delineated. In consequence, by considering that the magnitude of its v x q and v y o 
components range from km s _1 to the value given by this factor, all types of orbits are 
included. Thus, the initial conditions are just a point in the parameter space of any type of 
orbit of arbitrary eccentricity, size and orientation. 

The final step in the modeling consists of an iterative fitting of orbits to the observed 
spatial points (xrl , yiu) in the ridge-line of the intensity map (Fig. [2]) and the projected 
velocities (v z ). We evaluate the goodness of fit by means of an hybrid reduced \ 2 obtained 
by averaging the reduced x 2 °f eacn fitted parameter (xrl ,Vrl ,v z )- We weighted the x 2 
calculation by the square of the uncertainities - 1/4 of the spatial resolution for the x and 
y locations, and for v z the errors calculated during the extraction of the line properties 
as discussed in Section 13.2.11 If the model is a good approximation to the data, x 2 ~ 1- 
However, the hybrid nature of our x 2 an d the somewhat large uncertainities in the velocity 
values, could influence the scaling of the reduced x 2 making its exact value unimportant. In 
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this case, the best-fit will correspond to the orbit presenting the minimum \ 2 - 

The first attempts to model the kinematics showed that a single Keplerian orbit cannot 
reproduce the totality of v z vectors of the observed velocity field. Therefore, we decided to 
investigate two types of sets of initial conditions as suggested by the observed blue-shifted 
and red-shifted kinematical components: one for the northern region with xq = ± 5 pc, 
Do = 25 ± 5 pc and v z q = —25 ± 10 km s -1 ; and one for the southern part with x = —10 ± 5 
pc, i/q = —25 ± 5 pc and v z0 = 30 ± 10 km s -1 . The initial spatial coordinates (x ,y ) 
of the two sets were selected based on the morphology and kinematics. At first, any good 
fit to the data should follow the ridge-line. Therefore it is reasonable to have the north 
and south starting points located on it. Particularly for the northern region, the initial 
coordinates correspond to the peak of H 2 emission in Fig. [2J Furthermore, they are located 
in regions where the velocity uncertainties are low. This is once more particularly important 
for the northern region where the blue-shifted velocities associated morphologically with 
the northern tongue are well recognizable. Finally, a visual inspection of the velocity field 
suggests that the motions closer to the nucleus are probably perturbed by other physical 
processes and could possibly mislead the modeling. 

Approximately 5 x 10 4 Keplerian orbits of arbitrary eccentricity, size and orientation 
were modeled for each region and compared to the observations to determine a best fit. The 
orbits providing the best approximations to the data at these scales can be further compared 
with the large scale intensity map shown in Fig. [T] in order to study the connection between 
the gas in the central arcsecond and the circumnuclear environment. Particularly, we will 
be able to test the hypothesis of gas streaming from the circumnuclear ring along the linear 
structure. This hypothesis could motivate a different approach for the determination of the 
initial conditions (x , y , v z0 ). The implicit assumption here is that the gas observed in the 
central arcsecond (basically the southern and northern tongues) originates from the ring and 
thus the initial conditions are determined by the crossing points of the linear structure with 
the ring. However, we decided not to follow this approach because the crossing points of 
the linear structure with the ring are not clearly determined, especially for the southeastern 
region where several filaments and knots of gas are connecting the ring to the linear structure 
(Fig. [T]). In a ddition, the complex velocity field at these scales contains several kinematical 



components (jMuller Sanchez et al.l l2008bl ). This complicates greatly the determination of 
the initial v z vectors. Therefore, the approach based on the two sets of initial conditions 
mentioned in the last paragraph provides a reasonable more general methodology, which 
can test the hypothesis of gas streaming from the ring along the linear structure, and also 
determine the exact crossing points of the linear structure with the ring. 
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3.2.3. Gas is on highly elliptical/parabolic orbits 

Based on the reduced x 2 goodness of fit analysis described in Section |3.2.2[ we found for 
the northern and southern regions that among the complete set of tested orbits, only highly 
elliptical/parabolic orbits provide good approximations to the data points (xrl ,Vrl ,v z ). 
These orbits are contained in planes that are rotated through an axis coinciding with the 
x— axis (PA of the rotation axis is 90°) and have inclinations between 30 — 60°. These results 
confirm the presence of gas streamers in the nucleus of NGC1068. The initial conditions and 
other parameters of three representative orbits of the southern region are shown in Table [TJ 
The same information for three orbits of the northern region is presented in Table |2J The 
trajectories of the six orbits in the projected plane of the galaxy are shown in Fig. [7] plotted 
over the intensity map. As can be seen in the left panel of this Figure, the southern orbits 
are spatially associated with the southern tongue, and thus they are identified in Table [1] 
as orbits/streamers ST1, ST2 and ST3. The northern orbits are clearly associated with the 
northern tongue (right panel of Fig[7j), and therefore they are designated as orbits/streamers 
NT1, NT2 and NT3. 

The v z components of the six orbits are presented in Fig. [HJ For comparison, velocity 
curves of the data were extracted along the trajectories of the streamers. As these velocity 
curves are all similar, only the ones corresponding to orbits ST2 and NT2 are plotted in 
the left and right panels of Fig. [H] respectively. As can be seen in the two panels of this 
Figure, there is good agreement between the models and the data within the error bars. 
The velocity of the gas that is now lying in front of the nucleus corresponds to that of the 
southern streamer. This is demonstrated in both panels of Fig. [8] and in the left panel of 
Fig. E] which shows the trajectories of orbits ST2 and NT2 plotted over the velocity map. 
It can be seen in this Figure that the southern red-shifted component is connected to the 
nuclear kinematic component along the streamer trajectory. For the northern region, the 
projected velocities of the models fit well the data until r ~ 0.17". Within this radius 
the velocities of the streamer are not recognizable as they are being hidden by the nuclear 
redshifted kinematical component formed basically by the southern streamer and motions 
resulting from cloud-cloud collisions at few parsecs from the nucleus. 

We evaluate the goodness of fit to the ridge-line and the projected velocities by means 
of the reduced x 2 analysis discussed in Section 13.2.21 The results are shown in Tables [T] 
and [2J During the fitting process, it became clear that the somewhat large uncertainities in 
the velocity values influence the calculation of the reduced x 2 ■ m consequence, the best-fit 
corresponds to the orbits presenting the minimum x 2 ■ This was found for orbits ST2 and 
NT2 for the southern and northern regions, respectively. The conclusion arising from the 
orbital fits is that in the central arcsecond of NGC1068 there exist gas streamers flowing 
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along the southern and northern tongues in a plane which is rotated around an axis with a 
PA= 90° and has an inclination of 45°. 



3.2.4- Streamers flow inwards in the plane of the galaxy 



The orientation of the planes of the orbits discussed in Section 13.2.31 is so that the 
top half is on the far side and the bottom half towards us. However, the same fits to 
{pCRL , Url , v z ) are obtained for orbits located in planes that have the opposite orientation 
to the one presented in Tables [H and [2] (negative inclinations). In other words, based only 
on the kinematic modeling, there are two sets of orbits that we cannot distinguish: those 
in which the gas is heading towards the nucleus, i.e. the gas is streaming inwards (positive 
inclinations), and those in which the gas is flowing away from it, i.e. the gas is streaming 
outwards (negative inclinations). Luckily, our data combined with previous observations 
provide enough evidence to break this degeneracy. 

As demonstrated by several authors, the NE radio jet falls on the near side and its SW 
counterpart on the far side of the galaxy disk (e.g. Gallimore et al. 1994). The NE jet exhibits 
a bending point in its way out of the galactic plane caused by a shock interac t ion between 
the 



2002 



et and a dense molec ular cloud (jGallimore et al.lll996l ; iTecza et al.ll200ll ; I Cecil et al 



Poncelet et al.ll2007l ). The direct evidence for the molecular cloud obstructing the way 



of the jet in our SINFONI data (Fig. [5]) indicates that the orientation of the plane containing 
the cloud must be compatible with that of the galaxy plane, i.e. the top half is on the far side 
and the bottom half towards us. In consequence, the gas transport is towards the nucleus, 
i.e. the gas is streaming inwards. 

We can compare this scenario wi th previous observations of t his galaxy. According to 
the "NGC1068 Ringberg standards" (iBland-Hawthorn et al.l 119971 ). the galaxy disk has an 
inclination of i = 40 ± 3° and a kinematic major axis close to 90°. A similar configuration 
was found from the kinematics of the stars observed at r = 3" (~ 200 pc) with SINFONI 



(IDavies et al.l 120071 ; iMuller Sanchez et al.l l2008bl ) . These authors found that a rotating disk 
model with a PA= 84 ± 4° and i = 45 ± 4° reproduces the stellar kinematics well. At scales 
down to r = 0.5", the stellar velocity field is compatible with the same kinematic major axis 
but the inclination is not well constrained. This situation is similar with what we have found 
for the gas, although our modeling favours an inclination of 45°. 

The orientation is consistent with the strong evidence across the electromagnetic spec- 
trum that the NE radio jet, ionization cones and the extended x-ray emission fall on the 
near side, and the SW counterparts of all these fall on the far side, and therefore hidden by 
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the galaxy disk te land-Hawthorn et al.lll997l ). Fi nally, the plane orientatio n of our model is 
also consistent with the recent H 2 observations of iGalliano fc AUoinl (120021 ) who reproduced 
the observed changes in the line profiles across the 200 pc central region by a two-component 
kinematical model consisting of a rotating molecular disk inclined 65° with this orientation 
plus an outflow. 

Thus, the resulting physical model consists of gas streamers located in the plane of the 
galaxy, and falling inwards from either side of the nucleus. The streamer in all cases ap- 
proaches the centre (pericenters of a few pc), but as far as we can tell with the present spatial 
resolution, it does not actually reach the point source (see Tables [1] and [2]) . However, cloud- 
cloud collisions at few parsecs from the nucleus may provide a mechanism through which 
orbiting gas can lose angular momentum and remain in the inner region, falling eventually 
to the centre. Additionally, it can be seen in the right panel of Fig. [HI that the streamers 
constitute a u-dip in the whole field and thus are cold, indicating that the gravity dominat- 
ing the overall flow velocity is stronger than the internal dispersion. The tidal stretching of 
the material during infall produces the streamers now observed as a linear structure. These 
motions are strong evidence that we are seeing, on scales down to a few parsec, how gas is 
being driven toward the AGN in NGC1068, and hence how the AGN is being fuelled. 



3.3. Mass of the infalling molecular material 



We now discuss the molecular gas masses of the southern and northern tongues. The 
lack of robust direct probes of the gas mass complicates the analysis. To mitigate the 
uncertainties, we have used several independent methods. 



Our first estimate is based on the dynamical mass obtained by iDavies et ajj (120071 ) 
from the stellar kinematics in the central arcsecond of the galaxy. A dynamical mass of 
Mdyn = 1.3 x 10 8 M within r = 0.5" (35 pc ) is given by these authors. Assuming a 



conservative 10% gas fraction ( Hicks et al. 20081 ) and that most of the gas mass within this 
radius is distributed in equal amounts inside the southern and northern tongues (see Fig. [2]) , 
each of these components would have then a mass of ~ 6 x 10 6 M . This value can be 
considered as a lower limit due to the small adopted gas frac tion, which probably is higher 
in the nuclear regions of Seyfert 2 galaxies (IHicks et al.ll2008l ). 



On the other hand, assuming that the gas particles in each of the two tongues had 
initially a spheroidal geometry, we can esti mate a dynamical ma ss directly from the dispersion 
velocity from the relation Md yn = 5a 2 r/G ( iBender et al.lll992l ). Considering a spheroid with 
radius R = 0.1" (7 pc) based on the approximate FWHM of the southern tongue, and 
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(r = 70 i 15 km s" 1 as measured from the dispersion map along the trajectories of the 
streamers (Fig. [6]), we find a mass of 6.5 ±0.1 x 10 7 M Q for each concentration of molecular 
gas. The obtained value represents an upper limit on the gas mass since this calculation 
implicity assumes that the gas is gravitational bound and this may not be the case. 

A third independent method to estimate the mass of the molecular gas in the region is 
from the H? 1-0 S(T) luminosity and a n appropriate conversion factor between this quantity 



and M gas . iMiiller Sanchez et al.l (120061 ) found from large aperture (3" or more) measurements 



1-0S(1) 



/M t 



gas 



2.5 x 



of actively star-forming galaxies (including NGC1068) a ratio of L 
10~ 4 L & /Mq. Although this ratio was obtained for conditions which are probably not met 
in the nuclear region of the galaxy and depends on fraction of gas on cloud surfaces that 
can be heated, we apply it simply because there is evidence for vigourous star form ation 
close to the AGN jThatte et alill997l : bavies et ai1l2007i : IMiiller Sanchez et al.ll2008bh and 
it can help us to make at least an approximate estimate of the total molecular gas mass. 
From the total H2 1-0 S(l) luminosity of the northern and the southern tongue we estimate 
the total gas mass of each component to be 6 x 10 7 M Q and 4 x 10 7 M respectively. The 
uncertainties in this approach correspond to the uncertainties in the conversion factor which 
has a lcr statistical uncertainty from the 17 galaxies in their sample of a factor of 2. However, 
this factor has additionally a systematic uncer tainty inherent to the p robab le overestimated 
gas masses of the galaxies in the table from IMiiller Sanchez et al.l (120061 ) obtained using 
a standard CO-H2 conversion factor (overestimation of factors 2-5). This would lead to 
masses for the nortehrn and the southern tongue as little as 6 x 10 6 M and 4 x 10 6 M 
respectively. A s imilar conversion factor between cold and warm gas mass is proposed by 



Dale et all (12005[ ) for NGC1068 M^ rm /M£° ld = 1 x 10~ 6 . After converting the H 2 1-0 S(l) 



luminosities to warm gas masses (IDale et al.l 120051 ). we obtain total gas masses which are 
very consistent with the values calculated by means of the ratio £i-os(i) A^gas- 

The X-ray irr adiation of molec ular c louds by the centr al X-ray s ource in NGC1 68 ha s 
been discussed by iMalonev I Jl997h and iMatt et al.l Jl997l ). see also iNeufeld et all Jl994T ). 
Maloney et al. found that an attenuating column density of 10 22 cm -2 between the central 
engine and the H 2 emitting molecular clouds in the central r ~ 100 pc accounts for the 
observed H 2 intensity. As the gas density is a function of radius, with h igher densit i es the 
further in, it is probable that the column density is increased at our scales. IMatt et al.l (119971 ) 
have shown that the opaque material that obscures our direct view to the central engine of 
NGC1068 is Compton thick with an attenuating column density of at least Uh = 10 24 cm" 2 . 
However, this does not imply that the H 2 clouds have precisely, or even approximately, this 
value. By assuming this and 10 23 cm -2 as the upper and lower limits of the averaged column 
density in the central r = 0.4", we find that the southern/northern tongue has a gas mass 
ranging from 2 x 10 6 to 2 x 10 7 M Q . 
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Our last method for determining the gas mass es is based on the dus t emissivity at mid- 
IR wavelengths. Recent mid-IR observations from iTomono et al.l (120061 ) inferred a value of 



4.6 x 10 5 M m for the dust mass in the northe rn tongue. As the total gas mass is exp ected to 



be ~ 100 times greater than the dust mass JContini fc Continilliooa lDrainell2003h . a value 
of a 4.6 x 10 7 M can be inferred for the gas mass in this component. Assuming once more 
that the two gas concentrations have basically the same mass, the southern tongue will have 
also a gas mass of ~ 4.6 x 10 7 M . 

We can compare the values derived from the five different methods with previous esti- 
mates of the g as mass in the nuclea r region of NGC1068. From millimeter/sub-millimeter 
interferometry, iTacconi et al.l (11994] ) estimated the total mass contained within 1" of the 
nucleus as ~ 1.6 x 10 8 M M . More recent interferometric observations of the 12 CO(1-0) and 
12 CO(2-l) emission by ISchinnerer et al.l (120001 ) give a molecular mass of ~ 5 x 10 M Q con- 



taine d in the central 10 pc o f the galaxy. In addition, radiative transfer theoretical models 
from ISchartmann et al.l (120051 ) predict a dust mass of 8 x 10 4 M in the central 70 pc, which 
converted to gas mass results in a value of 8 x 10 6 M Q . 

All estimates are plausibly consistent with each other and with previous approximations 
suggesting that each molecular gas concentration (the southern or the northern tongue) has 
a mass within the range 6 x 10 6 — 6 x 10 7 M with a logarithmic mean of 2 x 10 7 M 
and a statistical uncertainty of a factor of 3. We will adopt this value for the mass of 
each component. This is consiste nt with a 25% gas f r action of the dynamical mass in the 



central arcsecond of the galaxy (ITacconi et a 



2000: Davies et al. 2007: Miiller Sanchez et al. 



19941: iThatte et all Il997i : ISchinnerer et al. 



2008bl ). Furthermore, it is con sistent with 



the molecular mass estimated from 12 CO(2-l) emission (ISchinnerer et al.ll2000l ) on slightly 
larger scales and an averaged column density in the central ±0.4" of 10 2 



i24 cm 2_ 



3.4. Mass accretion rate 

The morphology and kinematics of the gas are strong evidence that we are witnessing, 
on scales of a few parsecs, how gas is being driven toward the AGN in NGC1068, and hence 
how the nucleus is being fuelled. 

The mass accretion rate down to a few parsecs from the AGN can be estimated assuming 
that material falls into the nucleus through the linear structure. The infalling time scale, 
defined as the time a gas cloud takes to travel from the initial position to the pericenter, is 
obtained directly from the modeling and has a value of 1.3 Myr for the northern streamer 
(see Tabled]). This factor and the total gas mass yield a mass accretion rate at these scales of 
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~ 15 M yr _1 with a la uncertainty of a factor of 3. This value is an upper limit since not all 
the gas flowing towards the nucleus will actually stay there. A similar value is found for the 
southern streamer assuming that its mass corresponds t o the mass that curr ently is located 
in front of the nucleus. From millimeter interferometry, iTacconi et al.l (11994] ) estimated this 
influx to be a few M yr _1 from the total mass contained within 1" of the nucleus (~ 1.6 x 10 s 
M Q ) and radial velocities (~ 50 km s _1 ), which is of the same order as our data. 

We can compare this inflow rate on scales of ~ 10 pc to that onto the black hole itself. 
The mass accretion rate at scales down to one Schwarzschild radius Rs can be estimated 
from the mass-to-luminosity conversion efficiency o f a black hole L = n c 2 dM/ dt, where rj is 
the accretion efficiency with typical values 0.1 — 0.3 (jEardley fc Presslll975l ) and L is the bolo- 
metric luminosity of the AGN, which for NGC1068 is ~ 8 x 10 44 ergs s -1 (ITelesco fc Harper 
19801 ). This yields a mass accretion rate of 0.03 — 0.09 M Q yr -1 , indicating that the dM/dt 
is reduced » 1000 times from r = 1 pc to a few Rs- An analogous situation is observed at 
our own Gala ctic Center, for which m ass accretion rates of 10~ 3 "~ 4 M yr _1 at r = 1 pc 
are observed (IGenzel fc Towned 1 19871 ). and dM/dt ~ 10~ 7 '" -8 M yr _1 at a few Rs from 
the point source S grA* are estimated for a bolometric luminosi ty of the point source of a 
few 10 38 ergs s _1 (jOzernoy fc Genzell Il996l ; iMarrone et al.ll2007l ). The apparent dM/dt is 
a strong function of radius, with much reduced dM/dt closer to the nucleus. Qualitatively 
this is a natural consequence of inefficient angular momentum transport. While this may 
be a coincidence, we speculate that the gas accretion mechanisms are similar for active and 
non-acti ye supermassiye bla ck holes. Recent simulations of accretion of stella r winds on to 



SgrA* (ICuadra et al.l 120061 ) and supermassive black holes in Seyfert nuclei (jSchartmann 



20071 ) , reveal that the cold gas streamed down to scales approximately ten times smaller and 
settled into a very turbulent disk. This suggests that indeed the processes both, helping and 
hindering gas inflow, are the same for Seyfert galaxies and the Galactic Center. 



3.5. Obscuration by inflowing gas 



We now discuss the implications of the presence of molecular gas in front of the nucleus 
to the obscuration of the AGN in NGC1068. First, we calculate the gas mass surface density 
S of the southern tongue assuming once again a cloud radius of R = 0.1". For a gas mass 
of 2 x 10 7 M we obtain £ = 12 x 10 4 M pc -2 , yielding a column density of n H = 8 x 10 24 
cm' 2 , which is comparable with the values predicted for highly absorbe d objects (tin > 10 
cm" 2 ) and speci fically for Compton thick sources such as NGC1068 (Bas sani et a.. 



Matt et al.l 119971 ) and the column densities of clumpy torus models (INenkova et al 



l24 



1999 



2002 . 



This high column density suggests that the nuclear structure is optically thick in the near- 
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IR and hence this gas concentration can be associated with the obscuring material that is 
hiding the broad line region. In addition, as this value represents an average on scales of 
~ 0.2", the true column density of individual smaller clouds must be larger. The fact that 
a considerable amount of the 2 /im and 10 /zm emission from the innermost region (~ 1 — 2 
pc) and the maser disk radiation is not obscured despite the presence of large gas column 
densities, implies that the structure must be a clumpy medium. Line-of-sight attenuation 
of all these in this case would be a mere consequence of one or more intervening molecular 
clouds. 

Our observations therefore suggest that the southern tongue can be associated with 
the obscuring material that is hiding the nucleus but not in the classical picture of a ro- 
tating torus. This scenario is mainly ruled out by the kinematics which do not show any 
type of rotation near the AGN. However, there are several pieces of evidence that sup- 
port this association. First, th e size scale of the nuclear gas is remarkably similar to those 



Tomono et al. 2006; Poncelet et al 



of recent mid-IR observations ( Bock et al.l |2000| ; iTomono et al.1 l200ll ; iGalliano et al.l 12005 



2007) and those of static torus models, in particular the 



latest clumpy model of Honig et al. (12006 ). for which a size of 15x7 pc (diameter) is pre- 



dicted for the H 2 distribution in this galaxy (Honig personal communication) . There is also a 
similarity between the PA of the major axis of the core ~ 120°, which is consistent with that 



of th e line of maser spots (jGreenhill et al.l Il996l ) and the 300K dust emission ( jJaffe et al. 



20071 ). Furthermore, the estimated gas mass and column density fully agree with previous 
observations and torus models (see Sect. 13. 3j) . Hence we interpret this nuclear concentration 
of set of infalling clouds, that form the optically thick outerpart of an amorphous 

clumpy molecular /dusty structure. This large scale structure that we have observed will 
most probably enclose smaller clouds, qualitatively similar to a nested clouds scenario in 
which a distinct rotating molecular /dusty torus may or may not be present. In any case, 
based on the morphology and kinematics of the gas, we can state that if there is a rotating 
torus in NGC1068, its outer radius R out must be smaller than 7 pc and it is encircled by this 
amorphous molecular /dust obscuration. 



3.6. Connection of the streamers with the circumnuclear environment 

The remaining subject to investigate is the origin of the infalling material. A detailed 
study of the physical conditions of the NLR and the molecular gas at larger scales is crucial 
for the understanding of this phenomenon. We postpone our analy sis of these and other 



featur es in the central 3" x 3" of NGC1068 to a subsequent publication (jMuller Sanchez et al. 



2008bl ). Here, we provide a few remarks on the possible connection of the streamers with 
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the circumnuclear environment. 

In the central arcsecond, the orbital fit shows that highly elliptical/parabolic orbits 
in the plane of the galaxy reproduce the ridge and velocities of the gas quite well. A 
visual inspection of Fig. [1] shows that apparently the central two streamers are part of a 
linear structure connecting the NNE and SSW to the circumnuclear ring. We can test 
this hypothesis by comparing the complete trajectories of the orbits with the larger scale H 2 
emission map in Fig. [TJ A streamer is consistent with an origin in the ring if its apocenter lies 
inside the ring. If the apocenter of the orbit lies outside the ring, a ring origin of the streamer 
is uncertain. In order to have a reference of the width of the ring, we have superimposed 
on the intensity map several concentric thin rings of different sizes but having the same axis 
ratio and position angle as the observed ring in Fig. [TJ Based on width of brightest H 2 
emission around the ring, this overlay suggests that one can assign a characteristic thickness 
to the ring of ~ 0.6" (~ 42 pc). 

The results of this investigation for the six orbits described in Tables [TJ and [2] are shown 
in Fig. [HI In this Figure, we have plotted the complete trajectories of the orbits over the larger 
scale H 2 flux map presented in Fig. [TJ However, as the gas is streaming inwards, we observe 
only the part of each orbit which is transporting gas to the nucleus (the eastern part of each 
orbit). Based on the results from the x 2 goodness of fit analysis, we considered primarily 
orbits with 45° inclination. As can be seen in Fig. [9j orbit NT2 follows the northern part of 
the linear structure until it apparently merges with the circumnuclear ring. The rest of the 
streamer's trajectory is fully contained in the ring implying that the streamer is consistent 
with an origin in the circumnuclear ring. In the southern region, orbit ST2 follows the nearly 
linear extension of the southern tongue until r ~ 0.6". After this point, the linear structure 
is not distinguishable anymore from the several filaments and knots of gas which appear to 
be emanating from the ring. In this the apocenter of the orbit is contained in the 

ring (note that the width of the ring can be slightly larger), streamer ST2 is also consistent 
with an origin in the circumnuclear ring. 

The previous analysis has not only eludicated a plausible connection of the gas streamers 
with the circumnuclear environment, but also confirmed the results obtained from the \ 2 
goodness of fit analysis. On one hand, the large scale study suggests that orbit ST3 is not 
a good fit to the data at these scales as most of its trajectory is actually observed in the 
empty cavity of the emission map where very little gas is detected. This and the fact that 
the origin of orbit NT3 is uncertain imply that the gas is with high probability not contained 
in orbits with inclination angles close to 60°. On the other hand, if one asserts that the gas 
originates in the ring, orbits with 30° inclination (ST1 and NT1) provide a moderate fit to 
the data at these scales. Notice that due to the uncertainty in the thickness of the ring, 
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the origin of streamer NT1 could be associated with the ring. In any case, if the streamers 
originate in the ring, the large scale study clearly favours orbits with inclination angles close 
to 45°, fully consistent with the results from previous \ 2 analysis (Section 13.2.31) . 

Whether the circumnuclear ring is contained in the plane of the galaxy as the streamers 
or in another configuration is presently uncertain. Assuming a ring origin for the streamers, 
and that the c ircumnuclear ring is contained in the plane of th e galaxy as suggested by pre- 
vious authors (ISchinnerer et al.ll2000l : iGalliano fc Alloml 120021 ). then for material in the ring 
to approach the center, collisions must remove a significant fraction of the angular momen- 
tum. There is strong evidence that gas in the ring exhibits significant no n-circular motions 
jGalliano fe Alloinlbooa bavies et aDbood : Miiller Sanchez et alil2008bf ). This could result 
in collisions that lead to loss of sufficient angular momentum. Such non-rotating clouds 
would likely collide with orbiting material in the ring at a stationary point in space, the 
gas torn off the ring would originate at this stationary point, and so should follow a time- 
independent path. The tidal stretching of the material during infall produces the streamers 
now observed as a linear structure. However plausible this scenario may seem, it remains 
speculative at this point. A better understanding of the physical properties of the circum- 
nuclear ring is required. T hus, as mentioned above, we postpone any interpretation of the 
ring to a subsequent paper (jMiiller Sanchez et al.ll2008bl ). 



4. Conclusions 

We have presented in this paper high resolution SINFONI observations of the molecular 
gas in the nucleus of NGC1068. The distribution of the H 2 1-0 S(l) emission at a resolution 
of 0.075" has been resolved. Two bright regions connected by a linear structure extending up 
to 0.7" north of the nucleus along a PA of ~ —14° are distinguished: one lying right in front 
of the nucleus (the southern tongue) and another one 0.35" north of the center (the northern 
tongue). The northern tongue correlates with the mid-IR emission and its tip coincides with 
a knot of radio continuum emission providing direct evidence of the shock interface between 
the jet and a molecular cloud that has caused the jet direction to change. The main results 
of our analysis on the kinematics of these components are summarized as follows: 

• Dynamical modeling shows that material is streaming towards the nucleus. The in- 
falling gas is contained on elliptical/parabolic orbits whose orientation is consistent 
with that of the plane of the galaxy. We interpret this as strong evidence of how gas, 
on scales of a few parsecs, is fueling the AGN. 

• The gas transport is from ~ 70 pc to a few parsecs from the nucleus in the plane of 
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the galaxy. The modeling reveals the existence of two streamers: a northern streamer 
associated with the northern tongue that passes very close to the nucleus (pericenter 
of 5 pc), and a southern streamer which lies currently in front of the nucleus associated 
with the southern tongue, and has a pericenter of ~ 1 pc. 

• The mass inflow rate dM/dt is ~ 15 M yr -1 from scales of 30 pc to a few pc. This 
is about 1000 times that from a few pc to a few times the Schwarzschild radius R S - A 
similar change in the mass inflow rate with radius is observed in the Galactic Center. 
This may be a natural consequence of inefficient angular momentum transport. 

• The geometry, kinematics and high column density of the nuclear concentration of 
molecular gas (the southern tongue) can be explained by a tidally disrupted streamer 
consisting of a set of infalling clouds that form the optically thick outerpart of an 
amorphous clumpy molecular /dusty structure. 



The authors are grateful to the staff at the Paranal Observatory for their support during 
the observations, and to the entire SINFONI team at MPE and ESO. 
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Table 1: Orbital parameters of the southern gas streamers 



Orbit 


i a 


PA a 


xo b 


y h 


_ b 


Vxo c 


VyO ° 


Vzo c 


J/J, 


y P d 


z d 


x a 


Va 


z e 


t f 


x 2 


ST1 


30 


90 


-11 


-25 


-14 


6 


35 


20 


0.1 


0.42 


0.24 


-25 


-121 


-70 


0.9 


0.87 


ST2 


45 


90 


-11 


-25 


-25 


-1 


25 


25 


-0.01 


1.24 


1.24 


5 


-70 


-70 


1 


0.18 


ST3 


60 


90 


-11 


-25 


-42 


-1 


17 


30 


-0.3 


1.07 


1.85 


28 


-98 


-169 


1.1 


0.71 



inclination and position angles are given in (°) 
fc Initial spatial coordinates in (pc) 
c Initial velocity vectors in (km s _1 ) 
^Spatial coordinates of the pericenter in (pc) 
e Spatial coordinates of the apocentcr in (pc) 

^Infalling time scale in (Myr). This is defined as the time a gas cloud takes to travel from the initial coordinates 
to the pericenter 



Table 2: Orbital parameters of the northern gas streamers 



Orbit 


i & 


p A a 


x h 


y h 




Vxo c 


Vy C 


VzO C 


rr. d 


y P d 


z d 




Va 


z e 


t f 


x 2 


NT1 


30 


90 


0.0 


25 


14 


-25 


-35 


-20 


-5 


-3.7 


-2.1 


82 


61 


35 


1.3 


1.1 


NT2 


45 


90 


0.0 


25 


25 


-18 


-25 


-25 


-2.9 


-2.9 


-2.9 


56 


54 


54 


1.3 


0.51 


NT3 


60 


90 


0.0 


25 


42 


-12 


-20 


-35 


-2.4 


-1.8 


-3 








1.1 


2 



inclination and position angles are given in (°) 
fc Initial spatial coordinates in (pc) 
c Initial velocity vectors in (km s _1 ) 
^Spatial coordinates of the pericenter in (pc) 
e Spatial coordinates of the apocentcr in (pc) 

^Infalling time scale in (Myr). This is defined as the time a gas cloud takes to travel from the initial coordinates 
to the pericenter 
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Fig. 1. — SINFONI 2.12 /im H 2 1-0 S(l) flux map of the central region of NGC1068 at 
a reso lution of 0.1". A contour plot of the 12 CO(2-l) emission map from ISchinnerer et al. 
(120001 ) is overlaid on the H 2 flux map. The scale is given in arcsec (1" = 70 pc). The position 
of the nucleus as defined by the near-IR non-stellar peak is represented by a crossed circle. 
In addition to the previously known prominent gas concentrations ~ 1" east and west of the 
nucleus, the morphology reveals several new distinct structures. One of them is a complete 
ring centered ~ 0.6" SW of the nucleus. The other is a linear structure apparently emanating 
from the circum- nuclear ring and connecting to the AGN along a PA of ~ —14°. This new 
linear streamer is the focus of the discussion in this paper. The central square emphasize the 
location of this linear structure and denotes the field of the small scale observations shown 
in Figures 2-8. 
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Fig. 2. — Flux map of the H2 1-0 S(l) emission in the central ±0.4" of NGC1068 with a 
resolution of 0.075" (~ 5 pc). The peak of the non-stellar con tinuum is represented by a 
crossed circle. The image is binned using Voronoi tessellations (jCappellari fc Copinl 120031 ). 
The line properties could not be extracted in the very central region around the AGN due to 
the strong continuum emission. The bins where the line properties could not be extracted 
were masked out and correspond to the diagonal patterned regions. The open triangles show 
the projected trajectory of the northern concentration of gas (the northern tongue, Orbit 
NT2). The half-crosses show the past trajectory of the gas which is currently located in 
front of the AGN (the southern tongue, Orbit ST2). 
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Fig. 3. — Integrated spectrum of the nucleus of NGC1068 in the wavelength range between 
1.5 and 2.35 /im (H + K— bands). The individual spectra were added over a 0.3" x 0.1" 
rectangular aperture centered at the position of the nucleus. The dashed line represents the 
best fitting curve to the continumm emission which corresponds to a black body radiating at 
700 K. The dotted square emphazises the region between 2.05 and 2.20 /im, the wavelength 
range where the H2 1-0 S(l) emission line can be found. This region of the electromagnetic 
spectrum is plotted in the small panel top left. The continuum has been substracted in 
this panel. The H 2 1-0 S(l) emission line is clearly seen. Other relevant emission lines are 
indicated in the figure. 



-28- 



■2, 



-a -1 o 1 2 
RA off sat (arc sec) 





-0.2 



-0.4 

-0.4 -0.2 0.0 0.2 
RA offnt (arcMc) 



a.4 aoo I 



Fig. 4. — Contour plot of the 12.5 /im deconvolved image of iBock et al.l (120001 ) overlaid on 
the SINFONI H2 1-0 S(l) flux maps. The left panel shows the overlay in the central 4" x 4" 
(mid-IR image in white contours) and the right panel in the central 0.8" x 0.8" of NGC1068. 
The mid-IR centre is located at the position of the nucleus as defined by the near-IR peak. In 
the right panel the boundary of the northern tongue is denoted by dashed lines emphazising 
the good correlation between the two images. 
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Fig. 5. — Morphology of the H 2 1-0 S(l) emission with the 5 GHz radio continuum over- 
plotted in contours showing the position of the nucleus SI and the jet-cloud interaction at 
the position of component C. 
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Fig. 6. — Velocity (left) and dispersion (right) maps of the molecular gas extracted from the 
SINFONI datacube in the central ±0.4" of NGC1068. Velocities are measured with respect 
to the systemic velocity. The dispersion has been corrected for instrumental broadening. 
A crossed circle indicates in each cas e the peak of the contin uum emission. The maps 
are binned using Voronoi tesellations (ICappellari fc Copinll2003l ). The bins where the line 
properties could not be extracted were masked out and correspond to the white regions. The 
rejected pixels in both maps are those with a flux density lower than 10% of the peak of the 
central emission shown in Figure [2] and are shown in white in the right part of the fields. 
The open triangles show the projected trajectory of the northern concentration of gas (the 
northern tongue, Orbit NT2). The half-crosses show the past trajectory of the gas which is 
currently located in front of the AGN (the southern tongue, Orbit ST2). 
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Fig. 7. — Left: Trajectories of orbits ST1 (red curve), ST2 (blue curve) and ST3 (green 
curve) plotted over the intensity map in black-white colour scale. Right: Trajectories of 
orbits NT1 (red curve), NT2 (blue curve) and NT3 (green curve) plotted over the intensity 
map in black-white colour scale. The peak of the non-stellar continuum is represented by a 
crossed circle in both panels. The bins where the line properties could not be extracted were 
masked out and correspond to the diagonal patterned regions in both images. 
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Fig. 8. — Line-of-sight velocities (v z ) of the streamers. The left panel shows orbits ST1 (red 
dashed line), ST2 (blue dashed-dotted line), and ST3 (green dotted line). The right panel 
shows orbits NT1 (red dashed line), NT2 (blue dashed-dotted line) and NT3 (green dotted 
line). In the left panel, the black solid line corresponds to the velocity curve extracted from 
the data along the trajectory of orbit ST2 starting at r = 0.4" until r = 0.015" (1 pc). In the 
right panel, the black solid line corresponds to the velocity curve extracted from the data 
along the trajectory of orbit NT2 starting at r = 0.35" until r = 0.2" (14 pc). The jet/cloud 
interaction in the north takes place at r~ 0.25, and is indicated by a red circle in the right 
panel. 
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Fig. 9. — Complete trajectories of the six orbits described in Tables [T] and [2] plotted over 
the large scale H2 intensity map. Orbits with 30° inclination (ST1 and NT1) are represented 
by red curves. Orbits with 45° inclination (ST2 and NT2) are shown in orange. Orbits with 
60° inclination (ST3 and NT3) are represented by green curves. The location of the AGN as 
indicated by the peak of the non-stellar continuum is represented by a white crossed circle. 
The initial spatial coordinates (xq and yo) of the orbits are indicated by small white circles. 
The white inner and outer concentric ellipses delinate the thickness of the molecular ring. 
The central white ellipse is only plotted for reference. 



